function [f,g]=fun2(x); 
f=100*(x(2)-x(1)^2)^2+(1-x(1))^2; %原函数
df=[-400*x(1)*(x(2)-x(1)^2)-2*(1-x(1));200*(x(2)-x(1)^2)]; %一阶导数
d2f=[-400*x(2)+1200*x(1)^2+2,-400*x(1) %二阶导数
 -400*x(1),200]; 